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Abstract. - Two-component fermionic superfluids on a lattice with an external non-Abelian gauge 
field give access to a variety of topological phases in presence of a sufficiently large spin imbalance. 
We address here the important issue of superfluidity breakdown induced by spin imbalance by 
a self-consistent calculation of the pairing gap, showing which of the predicted phases will be 
experimentally accessible. We present the full topological phase diagram, and we analyze the 
connection between Chern numbers and the existence of topologically protected and non-protected 
edge modes. The Chern numbers are calculated via a very efficient and simple method. 



C^ 



H Introduction. — The Landau theory of symmetry 

' breaking constitutes the fundamental tool to character- 

yize the different phases of matter. The discovery of the 



of realizing all possible classes of topological insulators in 
a single experimental setup (3J|8]. Artificial gauge fields 
for neutral atoms may be realized via the aid of, as ex- 
amples, adiabatic Raman passage [9], adiabatic control of 



superpositions of degenerate dark states 10 11 , spatially 



Q Quantum Hall (QH) and Quantum Spin Hall (QSH) effects 

C^ showed however, that some physical systems undergo tran- 

^"sitions between distinguishable phases without breaking varying Raman coupling 12 
fvq any local symmetry [I] . These transitions have a topolog- 

J>ical character, in the sense that the corresponding "order cently Abelian fields have been successfully realized exper- 
C^^ parameter" is not a local quantity, but rather an integer 
^^ number describing the system as a whole. 



tions to auxiliary states in optical lattices 13 ■ 18 



or Raman-induced transi- 
Re- 



OO 

■^ A strong theoretical and experimental effort is now fo- 
[^ cusing over topological systems, since those are known to 
C^ contain anyons, quasi-particles which are neither bosons 
^^nor fermions. These excitations may have non-Abelian 
\ , statistics, and are intrinsically robust and protected 
!>• against decoherence and external perturbations. Non- 
L^ Abelian anyons have been proposed as the key ingredient 
''j_J for the realization of a topologically-protected quantum 
d computer [2] . Systems with topological phases include var- 
ious classes of topological insulators [s] , the v = 5/2 Frac- 
tional Quantum Hall (FQH) state, and two-dimensional 
(2D) polarized fermionic p-wave superfluids pip]. 



imentally 19-21 , and non-Abelian fields may be available 
soon. 



In this paper we study a two-component fermionic gas 
on a 2D square lattice, with an external non-Abelian 
gauge field effectively yielding a spin-orbit (SO) coupling 
in the gas. We take a vector potential of the form 
A — {aUy, Pax), where ai are the Pauli matrices, and a, /? 
are independently tunable parameters. This potential can 
be generated via the techniques proposed in [7 18 22 . In 
the absence of interactions, the excitations of this system 
include massless Dirac fermions, and the associated Dirac 
cones may be tuned to be asymmetric in the x and y direc- 
In presence of an additional Abelian field, the 



tions 



23 



ground state is the so called squeezed Laughlin state, and 



Ultracold gases in optical lattices are at the forefront of 
current research in the quest for challenging many-body 
states [6|. Particularly interesting appear the recent pro- 
posals for the realization of artificial gauge fields for neu- 
tral atoms "t], since these offer the appealing possibility 



the system exhibits anomalous Integer QH physics 23 



Interactions are fundamental in order to access both the 
superfluid regime and the FQH physics. We study here 
the case of attractive interactions, which at low temper- 
atures naturally yield BCS superfluid pairing. In pres- 
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ence of a strong external magnetic field the superfluid 
enters a topological phase characterized by anyonic ex- 
citations with non-Abclian statistics 24 25 . The corre- 



obtain the complete Hamiltonian: 



sponding Hamiltonian is unitary equivalent to the one of a 
p-wave superfluid 124 . Nonetheless, the imbalance needed 
to observe a phase transition towards a topological state 
is larger than the one predicted to yield pair breaking, the 



so called Chandrasekar-Clogston (CC) limit 26 



We investigate here this important and still unaddressed 
question: is a topological phase accessible at all in ultra- 
cold s-wave superfluids? As it is shown in the following, 
we find a positive answer to this question: the spin-orbit 
coupling makes a superfluid sufficiently stable to enter var- 
ious topological phases. We continue by presenting the full 
topological phase diagram, calculated as a function of the 
lattice filling and the strength of the spin-flipping cou- 
plings along the x and y directions. The phase diagram is 
characterized in terms of the corresponding Chern num- 
bers (CN), which we obtain via a peculiarly simple and 
efficient method. Topological phases may be efficiently 
discussed in terms of the symmetry classes introduced by 
Altland and Zirnabuer 
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In our system Time-Reversal 
and Spin-Rotation invariances are destroyed by the Zee- 
man and SO terms, and as a consequence our Hamiltonian 
belongs to most general symmetry class D. The periodic 
classification of topological insulators [3[[8] predicts that 
the topological phases of this Hamiltonian in 2D are in- 
dexed in terms of a Z number. As we will see, phases 
with CN=0, ±1, ±2 are accessible with ultracold fermionic 
gases in presence of the strong SO coupling provided by 
a synthetic gauge field. We study different classes of edge 
modes in this system, evidencing their strict relation with 
the Chern numbers. In analogy with the QSH effect, we 
also find that the topological protection of the edge modes 
is destroyed when pairs of states with same spin appear 
on the same edge. Topological protection of these modes 
is therefore characterized in terms of the Z2 index given 
by the parity of CN. 

Description of the system. — We consider a spin- 
imbalanced mixture of two fermionic species with inter- 
species attractive interaction. Particles are free to move in 
a 2D square lattice in presence of an external non-Abelian 
field. Introducing c\ 
may be written as 



(C;^, 



Cj, ), the hopping contribution 



n 



KIN+SO 



-tX! H (cl_^_-7i^i+eCi + h.C, 



The hopping terms are defined as: 



li^i+v = 



(1) 



(2) 



Using the identity e*-^"^^ = ±2x2 cos A + iaj sin A, it is 
easy to see that the spin-conserving [spin-fiipping] hopping 
terms in the x and y directions are proportional respec- 
tively to cos(a) and cos(/3) [sin(Q;) and sin(/3)]. Adding 
the chemical potential, Zeeman and interaction terms, we 



i 

- ^ [(Ai + h)cyi^ + {^l- h)clci^'\ . (3) 



Here F > is the attraction strength, the BCS s-wave 
pairing Ai = —V{ci^Ci^) has been taken real, and a con- 
stant term has been neglected. Due to the mean-field char- 
acter of the BCS approximation, it is important to realize 
that the pairing strength A; is not an external parameter 
which is tunable at will, but it is a quantity that has to be 
calculated self-consistently. This point will be analyzed in 
detail in the next Section. 

In presence of periodic boundary conditions, the pair- 
ing gap may be taken as constant throughout the system 
(Ai = A) and the Hamiltonian can be easily diagonal- 
ized 24]. We introduce ^j^. — (c|^ ^, c|^ , ,c_k,t,c_k,^) with 



c|^ = y-1/2 ^. e«kict, and obtain 



Hk = 



Ek - hcr^+gk- a- 



lACTy 

-ek + ha^ + gk • cr" 



(4) 



(5) 



where cr = (ax-,<yy)- The spin-conserving dispersion is 
given by 



ejc ~ — 2t(cos a cos fc^, -t- cos /3 cos k^ 



^J■, 



(6) 



while the spin-flipping term arising from SO coupling 
reads 



gj^ = 2< (sin /3 sin ky , — sin a sin k^ ) . 
The eigenvalues of this Hamiltonian are 



(7) 



Xi = ei + h' + \g^\' + A'±2^h^el + A^) + el\gu\^. (8) 

In presence of pairing (A 7^ 0), the spectrum is 
generically gapped, but the gap may close at the four 
distinguishable momenta which yield gko = 0, i.e., 
ko G {(0, 0), (0, tt), (tt, 0), (tt, tt)}. The four values of the 
spin-imbalance at which the gap closes are given by 



(9) 



Transitions between phases with different topological 
properties may only happen at gap-closing points. Since 
the balanced superfluid {h = 0) is known to be topolog- 
ically trivial, the latter condition states that to access a 
topological phase one has to consider imbalances satisfying 
the condition h > mink(,(ft.ko)- 
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Self-consistent calculation of the pairing gap A. 

— It has been recently pointed out 24, 25, 28 that topo- 
logically non-trivial phases should appear in fermionic s- 
wave superfluids in presence of both SO coupling and spin- 
imbalance. In these works, the pairing gap was introduced 
as an external parameter tunable at will. This assumption 
should be taken with care at large spin-imbalance, since 
it is well known that superfluidity breaks down when the 
Zeeman energy (proportional to h) becomes large as com- 
pared to Aq, the pairing gap a,t h = and T = 0. There, 
the free energy of the paired state becomes larger than the 
one of the normal (A — 0) state, and the system presents 
a flrst order phase transition from superfluid to normal. 
In absence of SO coupling (a = /3 = 0) and in the con- 
tinuum, the critical value is analytic and given by the so 
called Chandrasekar-Clogston (CC) limit, hcc = Ao/-\/2 
26 . It may be seen from eq. (|9| that /iko > hcc for ev- 
ery ko, a condition that excludes the possibility of having 
topological phases for fermionic s-wave superfluids in the 
absence of SO coupling. 

The self-consistent calculation of the pairing gap A pro- 
ceeds via the minimization of the free energy 



F 



7VA2 
V 



E 



k L 



fk 



ft' 

i=l 



In ( 1 -H e 



-Ak,./fcBT 



(10) 

Here Ak,i (i — 1, ... ,4) are the four solutions of eq. (IS]), 
and A^ is the number of lattice sites. In the absence of 
SO coupling (a = /3 = 0) and in the continuum limit, by 
solving eq. ( 10 ) we find that A drops abruptly to zero at 
h — hcc in agreement with the CC limit. The results of 
this calculation in presence of SO on a lattice are shown 
in fig. [l] where we plot the magnitude of A as a function 
of the spin imbalance h and the attraction strength V. As 
one cranks up the SO coupling, the sharp phase transition 
gets smoothened 29 , and we find that values of A > t 



are easily reachable in the regions where h > minko(/iko)- 
Largest gaps are obtained when the spin-conserving and 
spin-fiipping tunneling energies have comparable magni- 
tudes, i.e., for |tan(a)| ^ |tan(/3)| ^ 1. As we will see 
in the following section, all phases with |CN| < 2 may be 
reached at ft. > minko(/iko) for particular values of a,/3, 
and ^. In fig. [1] this imbalance yields respectively (from 
top to bottom) CN=l,-l,-2. The superfluid is instead al- 
ways unstable due to pairing breakdown at imbalances 
much larger than minko(/iko)- 

Topological phase diagram. Our Hamiltonian 

has built-in particle-hole symmetry, i.e. its eigenvalues 
come in ± pairs. The topological state of this system 
may then be characterized in terms of the bulk Chern 
number (CN) of the upper band of the spectrum. The 
latter is a Z topological invariant which may change only 
when the spectral gap between two bands closes. In ad- 
dition, the sum of CNs associated to an isolate set of 
bands is conserved before and after the band touching. 
We calculate the upper band's CN following the elegant 




Fig. 1: Pairing gap |A|/i as a function of imbalance h and 
attraction strength V. The lines represent respectively the 
CC limit hcc = Ao/\/2 (continuous), and the topological 
boundaries hoo (dotted), /iott (wide dashed), h-^o (dot-dashed) 
and h^^ (fine dashed), as defined in eq. H. The gauge field 
strengths and imbalance here are a = 7r/4, and (/3, fi) — 
(7r/4, -3t) (top), (37r/4, -3t) (center), and (7r/4, -^f/2) (bot- 
tom). The numbers inside the plots give the CNs of the most 
accessible phases. The Chern number vanishes when A = 0, 
since a gas in the normal phase is always in a trivial topological 
rState. 
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method proposed by Bellissard in ref. 30 . Let's assume 



that the gap closes at point ko as h crosses the critical 
value /iko ■ The Hamiltonian % has then two eigenvectors 
l^g ), I^Q ) with vanishing eigenvalue. In the neighbor- 
hood of (ko,/iko)5 we may approximate the Hamiltonian 
with an effective 2x2 matrix obtained by projecting % 
on Jt/iq ), IV^Q ). Using the Pauli matrices a = {(7x^ (Ty,(Tz), 
the effective Hamiltonian may be written as 



•Hcff (k, h) = ^(k, h)+a- /(k, h). 



(11) 



If the spectrum features a linear Dirac cone at (ko,/iko)i 
the Jacobian matrix J? has a non-zero determinant at this 
point. In this case, as h increases above /ikoj Hic change 
in Chern number of the upper band is given by the sign of 
the determinant itself. This quantity is called the Berry 
index of the gap-closing point. More generally, if at a 
given value h the spectrum features multiple Dirac cones 
centered in kg , kp , . . . , the change in Chern number at 
h — h is given by the sum of the Berry indices at all the 
touching points, 



ACN(/i) = ^sign{det[J^-(/i,k[,*0]}. (12) 

i 

From the latter formula, it is clear that each Dirac cone 
changes the Chern number of a band by ±1. If the deter- 
minant of the Jacobian at a given gap-closing point van- 
ishes, the Chern number of a band changes there instead 
by 0, or ±2 
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The latter case is however rather patho- 
logical. In the present context, it is verified only when 
either a or /? equal integer multiples of tt/2. 

In order to investigate the topological character of var- 
ious phases, we show in fig. [2] the curves corresponding to 
the four critical imbalances ft-ko defined in eq. ^. The 
curves are plotted as a function of /3 for fixed a — 7r/4, 
and for four increasing values of the chemical potential /z 
(from left to right). We consider only negative chemical 
potentials, since the Hamiltonian has built-in particle-hole 
symmetry, and reversing the sign of /i produces symmetric 
phases to the ones found here. 

At small imbalances an s-wave superfluid is topologi- 
cally trivial for any fx, and CN=0. Each time the im- 
balance h crosses one of the critical curves /ikoj the gap 
between the bands closes at the corresponding momen- 
tum ko, and a topological phase transition occurs. At 
low lattice fillings, we find phases characterized by Chern 
numbers CN=-1,0,1. 

For higher lattice filling there appear phases with Chern 
numbers as large as ±2. As we will see in the following 
section, where we will discuss edge states, the phases with 
CN=-1,1 are topologically protected, while the phases 
with CN=-2,0,2 are not protected. The topological pro- 
tection of the system is therefore characterized by the Z2 
number given by mod2(CN), in analogy with the theory of 
the Quantum Spin Hall effect 31 . At half-filling, fi = 0, 



As we have seen in the previous Section, not all regions 
of the phase diagram plotted in fig. [2] are accessible due 
to pair-breaking. Nonetheless, phases with equal CN are 
topologically equivalent, and we have shown that phases 
with all CNs (-2,-1,0,1,2) may be realized with sufficiently 
large A(~ t) when h > miuko ( | /iko I ) ■ 

Spectrum on a cylinder. — A characteristic feature 
of topological phases is the presence of intra-gap modes 
localized on open boundaries (edges) of the system. Ul- 
tracold gases with sharp edges may be realized along the 
lines proposed in ref. 
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It should also be remembered 
that the core of a vortex is topologically equivalent to an 
edge which has been closed on itself. The topologically- 
protected edge modes we describe here are therefore equiv- 
alent to the Majorana fermions known to exist, e.g., in the 
core of vortices of p-wave superfluids. 

To investigate this aspect, we calculate the spectrum of 
the Hamiltonian (p]) on a cylinder, i.e., respectively with 
open and periodic boundary conditions along x and y. 
Due to the broken translational symmetry along the x- 
direction, we diagonalize the Hamiltonian by performing 
a Fourier transformation along the y-direction only. 

In figs. ISlandHlwe show the spectrum of 'Hx,k as a func- 
tion of ky for various combinations of fj,, /3, and h. The 
spectrum on a torus would be generally gapped, but when 
open boundary conditions are taken into account isolated 
states appear inside the gap. The eigenvectors correspond- 
ing to these isolated eigenstates are linear superpositions 
of t and 4, spin components, exponentially localized on ei- 
ther of the boundaries, as may be seen in fig. [51 These are 
therefore termed "edge states". Edge modes in this con- 



the phase diagram is composed of regions with CN=-2,0,2, 
and is therefore not topological for any a, (3, and h. 



text have first been discussed in ref. 24 for the specific 
case a = j3 = 7r/4. 

Given the refiection symmetry of the spectrum with re- 
spect to the sign of /cj,, these states always appear in ±fc 
pairs. Edge excitations may cross at zero energy with 
a linear dispersion ai ky — Q or tt, similarly to lattice 
Dirac fermions. Their origin is topological, and these zero- 
energy crossings may appear or disappear only at the criti- 
cal points h — /iko where the bulk spectra of the upper and 
lower bands touch in zero, see e.g. fig.lsjil). Since the fully- 
balanced {h — 0) and the fully-polarized {h ^ t) states 
are not topological and do not support any edge modes, 
on a lattice we find that edge states are present only for 
imbalances satisfying minko(/iko) < h < maxko(/iko)- 

When a single pair of states is present, see e.g. figs. |3|i 
and[3p, the two counter-propagating edge modes with ±fc 
are localized on opposite boundaries. Since on a given 
boundary there is no state available for backward spin- 
conserving scattering, these states are topologically pro- 
tected. These conditions are realized when mod2(CN)=l. 

When instead the parity of CN is 0, as depicted in 
fig. |4j the system contains either zero or two pairs of 
edge states. At sufficiently low-energy, the eigenvectors 
at a given energy are characterized by symmetric quasi- 
momenta ky — ±A:, ±(7r — k). As may be seen in fig. ^ 
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Fig. 2: Topological phase diagram of the system as a function of /3, the gauge field strength in the y-direction. Here we have 
set a — 7r/4 and A = i. From left to right, the chemical potential is ^/{\/2t) = —2 (a), — 1.25 (b),— 0.5 (c), (d). The lines 
depict the critical imbalances /loo (dotted), /iott (wide dashed), /i,rO (dot-dashed), ft,r7r (fine dashed). The phases are indexed by 
the Chern numbers calculated from eq. (12 I: CN=0 (white), |CN| = 1 (light gray), |CN|=2 (dark gray), and the meshed phases 
are the ones with negative CN. The labeled dots in the first and third plots mark the imbalances plotted in figs. |3] and [4] 
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Fig. 3: Spectrum on a cylinder, plotted as a function of ky. In each figure, the gauge field strength j3 and the imbalance h are 
given by the labeled dots in fig. [2k. The lines mark the edge states, and the thick dots mark the position of the gap-closing 
points. The CN of the upper band is 1 (-1) in the left (right) image, while the central image depicts a topological critical point. 
The chemical potential is /i/(\/2i) ~ —2, and the energy on the vertical axis is in units of t. 
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Fig. 4: Same as fig.lsl but here the chemical potential is fi/{\/2t) = —0.5, while the gauge field strength /3 and the imbalance h 
correspond to the labeled dots in fig. |2|;. From left to right, the CN of the upper band is respectively -2,-1, and 0. The symbols 
+ and * identify the eigenstates plotted in fig. [5] 
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Fig. 5: Edge modes on a cylinder: continuous (dashed) lines 
depict the amplitude of the spin f (4-) component of the wave- 
function with a, /3, fi and h as in fig.Hfl). The momenta of the 
edge states are given by the symbols in fig.BTI). Here two coun- 
terpropagating modes with identical spin state are localized on 
the same edge. In analogy with the QSH effect, these states 
are not topologically-protected, since a weak perturbation may 
induce backscattering from one mode to the other. 



the counter-propagating edge states with ky = k and 
ky = —n + k are both localized on the same edge. When 
niod2(CN)=0, we see that either edge modes can perform 
spin-conserving backscattering, or that there are no edge 
modes at all. A phase with even parity of the CN does 
not contain therefore any topologically-protected states. 

It is interesting to note that fig. |4lni) depicts a physical 
state with two pairs of edge states but vanishing CN. To 
the best of our knowledge, this case has not been discussed 
so far in the context of cold atomic gases. 

Conclusions. — We have studied here interacting ul- 
tracold fermions in presence of non-Abelian gauge fields. 
We have demonstrated that spin-imbalance can lead to a 
variety of intriguing topological phases. We characterized 
these phases by the associated Chern numbers, we pre- 
sented the full phase diagram, and we discussed the link 
existing between edge states and Chern numbers. Fur- 
thermore, we have shown via a self-consistent calculation 
which of the phases can be realized experimentally, and 
which are non-physical due to pairing breakdown. In light 
of the recent ground-breaking experimental realizations of 
synthetic Abelian 
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fields in ultracold gases, and 
given the prospective to realize non-Abelian fields soon, 
the fascinating physics discussed in this manuscript should 
be readily accessible in the very close future. 
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